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Abstract 

We study the short-range nucleon-nucleon interaction in a nonrelativistic 
chiral constituent quark model by diagonalizing a Hamiltonian containing 
a linear confinement and a Goldstone boson exchange interaction between 
quarks. A finite six-quark basis obtained from single particle cluster model 
states was previously used. Here we show that the configurations which ap- 
pear naturally through the use of molecular orbitals, instead of cluster model 
states, are much more efficient in lowering the six-quark energy. 

I. INTRODUCTION 

Constituent quark models have been applied to the study of the nucleon-nucleon inter- 
action. In a category of such models the Hamiltonian contains a kinetic term, a confinement 
term and an effective one-gluon exchange (OGE) term. These models explain the short-range 
repulsion in the NN systems as due to the colour-magnetic part of the OGE interaction 
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combined with quark interchanges between the 3q clusters. Nevertheless, an effective meson- 
exchange potential, introduced through the coupling of mesons to 3q cluster collectively, was 
required in order to reproduce the intermediate- and long-range attraction (for a review see 
for example |j]-f§). 

Another category are the hybrid models [[§-[§]. There, in addition to the OGE interaction, 
the quarks belonging to different 3q clusters interact via pseudoscalar and scalar meson 
exchange. In these models the short-range repulsion in the iViV system is still attributed 
to the OGE interaction between the constituent quarks. The medium- and long- range 
attraction are due to meson exchange, as expected. 

In a recent exploratory work [|7|, by using the Born-Oppenheimer approximation, we 
calculated an effective NN interaction at zero separation distance, within the constituent 
quark model |8|-[10|. In this model the quarks interact via Goldstone boson exchange (GBE) 



instead of OGE of conventional models, and the hyperfine splitting in hadrons is obtained 
from the short-range part of the GBE interaction. An important merit of the GBE model is 
that it reproduces the correct order of positive and negative parity states in both nonstrange 
H and strange baryons |TD| in contrast to any OGE model. In Ref. [J7| we showed that the 



same short-range part of the GBE interaction, also induces a short-range repulsion in the 
NN system. Moreover, the long and middle range attraction of the NN potential will 
automatically appear due to the presence of a Yukawa potential tail in the qq interaction 
and due to 2ti (or sigma) exchanges. 

In |7j the height of the repulsive core was about 800 MeV for the 3 Si channel and 
1300 MeV for the 1 S'o channel. Such a result has been obtained from diagonalizing the 
Hamiltonian of Ref. || in a six-quark cluster model basis built from harmonic oscillator states 
containing up to two quanta of excitation. The six-quark states have orbital symmetries [6]o 
and [42]o, so that they contain configurations of type s 6 , s 4 p 2 and s 5 2s, with the centre of 
mass motion removed. In the flavour-spin space only the symmetries [33], [51] and [411] were 
retained. As shown in |7| they produce the most important five basis states allowed by the 
Pauli principle. Due to the specific flavour-spin structure of the GBE interaction, we found 



that the state |sV [42] [51]fs) was highly dominant at zero-separation between nucleons. 
The symmetry structure of this state implies the existence of a node in the nucleon-nucleon 
.S-wave relative motion wave function at short distances. This nodal structure will induce an 
additional effective repulsion in dynamical calculations based, for example, on the resonating 
group method. 

A central issue of the NN problem is the construction of an adequate six-quark basis 
states. In principle the choice of basis is arbitrary if a sufficiently large basis is considered 
in the Hamiltonian diagonalization. But, as in practice one considers a finite set, its choice 



is very important. Ref. |TJ] advocated the use of molecular-type single particle orbitals 
instead of cluster model-type states. These orbitals have the proper axially and reflectionally 
symmetries and can be constructed from appropriate combinations of two-centre Gaussians. 
At zero-separation the six-quark states obtained from such orbitals contain certain p n s 6 ~ n 
components which are missing in the cluster model basis. In Ref. |12] it has been shown that 
for an OGE model used in the calculations of the NN potential they lead to a substantial 
lowering of the lowest eigenstate, used in the calculation of the NN potential. The molecular 
orbitals have also the advantage of forming an orthogonal and complete basis while the 
cluster model (two-centre) states are not orthogonal and are overcomplete. 

Due to the predominance ( 93 % ) of only one component, namely |s 4 p 2 [42]o[51]i?5), in 
the ground state wave function obtained in a cluster model basis flU the GBE model is a 
more chalenging case to test the efficiency of a molecular orbital basis than the OGE model, 
where there is some mixture of states (see e.g. [f[|,|l2"|l). Here we show that by using molecular 
orbitals the height of the repulsion reduces by about 22 % and 25 % in the 3 Si and 1 S 
channels respectively. 

The paper is organized as follows. In Sec. 2 we briefly recall the procedure of constructing 
six-quark states from molecular orbital single particle states. In Sec. 3 we describe the GBE 
Hamiltonian || . In Sec. 4 we present our results for zero-separation NN interaction derived 
in the Born-Oppenheimer approximation for the IS = (10) and (01) sectors. The last section 
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is devoted to summary and conclusions. 



II. SIX-QUARK STATES FROM MOLECULAR ORBITALS 



Here we follow closely Ref. |TT| where the use of molecular orbitals in the construction of 
six-quark states was originally proposed, instead of commonly used cluster model states. Let 
us denote by Z the separation coordinate between the centres of the two clusters. At finite 
Z, in the simplest cluster model basis, each of the six quarks is decribed by an orbital wave 
function represented by a Gaussian centered either at Z/2 or —Z/2. These nonorthogonal 
states are denoted by R (right) and L (left) respectively 



i?(f) = *(f- Z/2), L(r) = ^{r + Z/2). 



(1) 



Alternatively, in a molecular basis we consider the two lowest states, a which is even and 
7i which is odd. These could be either the solutions of a static, axially and reflectionally 



symmetric independent particle model Hamiltonian (see for example [13]) or, as for the 
present purpose, can be constructed from R and L states. 

First we introduce pseudo-right and pseudo-left states r and I starting from the molecular 
orbitals a and tt as 



2~ 1/2 (a ± tt) for all Z, 



(2) 



where 



< r \r >=< l\l >= 1, < r\l >= 0. 



(3) 



On the other hand, starting from the cluster model states, one can construct good parity, 
orthonormal states for all Z by setting 



a 



71 



[2{l± < R\L >)}- 1 l 2 (R± L), 



(4) 



which, introduced in (2) gives 
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R + L 



± 



R-L 



1+ < R\L >yi 2 (1- < R\L >)V2 



(5) 



At Z — > one has er — > s and 7r — > p (with m = 0, ±1) , so that 



2 1/2 (*±P), 



(6) 



and at Z — > oo one has r R and Z — > L. 

From (r, Z) as well as from (er, 7r) orbitals one can construct six-quark states of required 
permutation symmetry. For the Sq symmetries relevant for the NN problem the transfor- 
mations between six-quark states expressed in terms of (r, I) and (er, tt) states are given in 



Table I of Ref . |TT[ . This table shows that in the limit Z — > six-quark states obtained from 
molecular orbitals contain configurations of type s n p 6 ~ n with n = 0,1,..., 6. For example 
the [6]o state contains s 6 , s 6 p 4 , s 2 p A and p 6 configurations and the [42] o state associated to 
the S'-channel contains s 4 p 2 and s 2 p 4 configurations. This is in contrast to the cluster model 
basis where [6]q contains only s 6 and [42] o only s 4 p 2 configurations []I4 ]. This suggests that 
the six-quark basis states constructed from molecular orbitals form a richer basis without 
introducing more single particle states. Here we examine its role in lowering the ground state 
energy of a six-quark system described by the Hamiltonian introduced in the next section. 



Using Table I of Ref. |0] we find that the six-quark basis states needed for the Si or 



1 S'o channels are: 

|33[6] [33] FS ) 

|33[42] G [33W 

|33[42] [51W 
|33[42] D [411] F5 ) 

42+[6] [33] FS 



1 

4 



V3 (s 4 p 2 



2 4 

s p 



[6] G [33] 



4 2 

S p 



4 2 

s p 



4 2 

s p 



2 4 

S p 



2 4 

S p 



2 4 



(s 6 - p' 

" ' n [42] [33] FS/ , 

[42]o[5lU). 
[42] [411] F5 ), 

^(s 6 + p 6 ) - (sV + sV)] [6] G [33] 



FS / ' 



F5 / ' 



(7) 
(8) 
(9) 
(10) 

(11) 
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42+[42] [33] FS 

42+[42] [51] FS 

42 + [42] [411] FS 
51+[6] [33] FS 



'1 



sY + sY [42] [33] FS/ . 

[42] [51] F5 ), 



4 2, 2 4 

s p + s p 



1 

4 



sY + sY\ [42] [411] FS 
V3 (s 6 - p e ) + V5 (sY - sY)] [6] G [33] 



FS i 



(12) 

(13) 

(14) 
(15) 



where the notation 33 and mn + in the Ihs of each equality above means r 3 £ 3 and r m £ n +r n £ m 
as in Ref. (see also discussion below). Each wave function contains an orbital part (O) 
and a flavour-spin part (FS) which combined with the colour singlet [222] q state gives rise 
to a totally antisymmetric state. We restricted the flavour-spin states to [33]^, [51]fs and 
[411] F 5 according to the discussion given in Sec. II of Ref. where the most important 
states have been selected by using a schematic version of the Hamiltonian introduced in the 
next section. 

In a cluster model, the most important basis states built from s and p harmonic oscillator 

states are 



S 6 [6] G [33] 



FS / ' 



*V[42] [33] 



*V[42] [51] 



OV-' J -\FS ' 



S y[42] [411] 



IFS 



(16) 



(17) 



(18) 



(19) 



These are the first four states given by Eq. (8) of Ref. 0. The fifth one, containing the 
configuraiton s 5 2s is not considered here. Its role in lowering the ground state energy by a 
few MeV proved to be negligible. Besides being poorer in s n p 6 ~ n configurations, as explained 
above, the number of basis states is smaller in the cluster model although we deal with the 
same [f]o and [/]fs symmetries and the same harmonic oscillator states s and p in both 



cases. This is due to the existence of three-quark clusters only in the cluster model states, 
while the molecular basis also allows configurations with five quarks to the left and one to 
the right, or vice versa, or four quarks to the left and two to the right or vice versa. At large 
separations these states act as "hidden colour" states but at zero separation they bring a 
significant contribution, as we shall see below. 

The matrix elements of the Hamiltonian (22) are calculated in the basis (7-15) by us- 



ing the fractional parentage technique described in Refs. ||1^J15[| and also applied in Ref. 



. A programme based on Mathematica |L6| has been created for this purpose. In this 



way every six-body matrix element reduces to a linear combination of two-body matrix el- 
ements of either symmetric or antisymmetric states for which Eqs. (3.3) of Ref. || can be 
used to integrate in the spin-flavour space. Then the linear combinations contain orbital 
two-body matrix elements of the type (ss \Vy\ss), (ss \Vy\pp), {sp\Vy\sp), (sp\Vly\ps) and 
(PP \Vj I PP)l = o wnere 7 — 11 , V 01 v'i see Eq. (25). Here we study the case Z = for which 
the following harmonic oscillator states are used 

| s > =7r- 3 / 4 /r 3/2 exp(-r72/? 2 ), (20) 
\p > = 8 l ^3-^ 2 n- l ^p- 5 / 2 r exp (-r 2 /2(3 2 ) Y lm . (21) 

In this basis the orbital two-body matrix elements of the linear confinement V con f = Cr 
potential (23) are calculated analytically (see Appendix D of Ref. 0). 

III. HAMILTONIAN 

The GBE Hamiltonian considered below has the form |J : 

H = J2m + j:§r- W^T + £ v °°*^ + £ v x(ra) , (22) 

i i ^'"'i z l~<i rn i i<j i<j 

with the linear confining interaction : 

V mn i(r i /) = -^-X c j Cr ij , (23) 
and the spin-spin component of the GBE interaction in its SUp(3) form : 



Vydnj) = { £ K(r«)Af Af 



,F=1 



+ £ ^(ry)Af Af + ^(r^OAfAf + fyfaW) a, • a„ (24) 

with A = \f2/3 1, where 1 is the 3x3 unit matrix. The interaction (24) contains 7 = tt, K, 77 
and 77' meson-exchange terms and the form of V 1 (r^) is given as the sum of two distinct 
contributions : a Yukawa-type potential containing the mass of the exchanged meson and a 
short-range contribution of opposite sign, the role of which is crucial in baryon spectroscopy. 
For a given meson 7, the exchange potential is 

= ll^-C"' - ''"K^ - ^eM-*\ - r f)}. (25) 
For a system of u and d quarks only, as it is the case here, the i^-exchange does not 
contribute. In the calculations below we use the parameters of Refs. M. These are : 



22 2 
h± = 9 JR = o.67, 9 4± = 1.206, 

All 4:71 4:71 



r = 0.43 fm, a = 2.91 /m~\ C = 0.474 fm' 2 , m u4 = 340 MeV, (26) 
/V = 139 MeV, n v = 547 MeV, ^ = 958 MeV. 
In principle it would be better to use a parametrization of the GBE interaction as given in |17| 



based on a semirelativistic Hamiltonian. However, in applying the quark cluster approach 
to two-baryon systems we are restricted to use a nonrelativistic kinematics and an s 3 wave 
function for the ground state baryon. With an s 3 variational solution the nonrelativistic 



Hamiltonian introduced above works generally well ||18|| . In particular, for the nucleon, the 
quantity < N\H\N > reaches its minimum at 969.6 MeV which is only about 30 MeV 
above the nucleon mass obtained in the dynamical 3-body calculations of Ref. |J. There 
the shifted Gaussian of Eq. (25) results from a pure phenomenological fit. 



IV. RESULTS 



We diagonalize the Hamiltonian (22) in the six-quark basis (7-15) and calculate the ./V7V" 
interaction potential in the Born-Oppenheimer approximation 
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V NN (Z) = (H) z - <#>oo , (27) 

where (H) z is the lowest expectation value obtained from the diagonalization at a given Z 
and {H)oo = is the energy (mass) of two well separated nucleons. Here we study the case 
Z = 0, relevant for short separation distances between the nucleons. In Tables I and II we 
present our results for IS = (01) and (10) respectively, obtained from the diagonalization of 
H. From the diagonal matrix elements Ha as well as from the eigenvalues, the quantity 2m n 
= 1939 MeV has been subtracted according to (27). Here mjy is the nucleon mass calculated 
also variationally, with an s 3 configuration, as mentioned at the end of the previous section. 



This value is obtained for a harmonic oscillator parameter j3 = 0.437 fm [19]. For sake of 
comparison with Ref. we take same value of (3 for the six-quark system as well. 

In both IS=(01) and (10) cases the effect of using molecular orbitals is rather remarkable 
in lowering the ground state energy as compared to the cluster model value obtained in the 
four dimensional basis (16)-(19). Accordingly, the height of the repulsive core in the £3 
channel is reduced from 915 MeV in the cluster model basis (see Appendix) to 718 MeV 
in the molecular orbital basis. In the 1 Sq channel the reduction is from 1453 MeV to 1083 
MeV. Thus the molecular orbital basis is much better, inasmuch as the same two single 
particle states, s and p, are used in both bases. 

The previous study [0, performed in a cluster model basis indicated that the dominant 
configuration is associated to the symmetry [42] q [51] fs- It is the case here too and one can 
see from Tables I and II that the diagonal matrix element Ha of the state |42 + [42]o[51] F s > 
is far the lowest one, so that this state is much more favoured than |33[42]o[51]i?s > . As 
explained above, such a state represents a configuration with two quarks on the left and four 
on the right around the symmetry centre. At Z — > oo its energy becomes infinite i.e. this 
state behaves as a hidden colour state (see e.g. Ref. |TJJ]) and it decouples from the ground 
state. But at Z = it is the dominant component of the lowest state with a probability of 
87 % for IS = (01) and 93 % for IS = (10). The next important state is |33[42] [51] F5 > 
with a probability of 10 % for IS = (01) and 4 % for IS = (10). The presence of this 
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state will become more and more important with increasing Z . Asymptotically this state 
corresponds to a cluster model state with three quarks on the left and three on the right of 
the symmetry centre . 

To have a better understanding of the lowering of the six-quark energy we present in Ta- 
bles III and IV the separate contribution of the kinetic energy KE, of the confinement V con f 
and of the GBE interaction V x to the dominant state in the cluster model |s 4 p 2 [42]o[51]fs) 
result and the dominant state in the molecular basis case respectively. Table III corresponds 
to the 3 Si channel and Table IV to the 1 S'o channel. We can see that V con f does not change 
much in passing from the cluster model to the molecular orbital basis. The kinetic energy 
KE is higher in the molecular orbital basis which is natural because the s 2 p 4 and p 6 con- 
figurations contribute with higher energies than s 6 and s 4 p 2 . Contrary, the contribution of 
the GBE interaction V x is lowered by several hundreds of MeV in both channels, so that 
E = KE + V con f + V x is substantially lowered in the molecular orbital basis. This shows 
that the GBE interaction is more effective in the molecular orbital basis than in the cluster 
model basis. Note that E differs from the value of the diagonal matrix elements of Tables I 
and II by the additional quantity 6m — 2m^r, where m = m u = m^. 

The practically identical confinement energy in both bases shows that the amount of 
Van der Waals forces, as discussed in 0, remains the same. However, the soft attraction 
brought in by the Van der Waals forces does not play an important role at short distances 
and it should be removed in further studies at intermediate distances. 

For both IS = (01) and (10) sectors we also searched for the minimum of (H)z=o as a 
function of the oscillator parameter (3. For IS = (01) the minimum of 572 MeV has been 
reached at /3 = 0.547 fm. For IS = (10) the minimum of 715 MeV was obtained at (3 = 
0.608 fm. These values are larger than the value of (3 = 0.437 fm associated to the nucleon, 
which is quite natural because a six-quark system at equilibrium is a more extended object. 
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V. SUMMARY AND CONCLUSIONS 



We have calculated the NN interaction potential at zero separation distance between 
nucleons by treating NN cLS cL SIX - quark system in a constituent quark model where the 
quarks interact via Goldstone boson (pseudoscalar meson) exchange. The orbital part of 
the six-quark states was constructed from molecular orbitals instead of the commonly used 
cluster model single particle states. The molecular orbitals posses the proper axially and 
reflectionally symmetries and are thus physically more adequate than the cluster model 
states. Due to their orthogonality property they are also technically more convenient. Here 
we constructed molecular orbitals from harmonic oscillator s and p states. Such molecular 
orbitals are a very good approximation |20] to the exact eigenstates of a "two-centre" oscil- 
lator, frequently used in nuclear physics in the study of the nucleus-nucleus potential. The 
problem of calculating an NN potential is similar in many ways. 

We have shown that the upper bound of the ground state energy, and hence the height 
of the repulsive core in the NN potential, is lowered by about 200 MeV in the 3 Si channel 
and by about 400 MeV in the x So channel. Hence using molecular orbitals is more efficient 
than working with a cluster model basis. A repulsive core of several hundred MeV is still 
present in both channels. Due to the specific flavour-spin symmetry of the GBE interaction 
the molecular type component |42 + [42]o[51]i?5 > becomes dominant at short range which 
implies that the NN relative motion S-wave function has a node at short distance due to 
the presence of the configurations s 4 p 2 and s 2 p 4 . The dominance of the [51]i?s symmetry 
will reinforce the repulsion in dynamical calculations. In fact, it has been shown |l| that 
the phase shift calculated within the resonating group method with a pure [51] ps state 
shows a behaviour typical for potentials with a repulsive core. In OGE models this effect 
is absent because none of the [42] o states is dominant (see e.g. Jl2|) . Note also that the 
configurations s 2 p or p G introduced through the molecular orbitals might have an influence 



on the momentum distribution of the NN system as was discussed, for example, in |2T 
within the chromodielectric model. 
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The following step will be to calculate the NN potential at Z ^ 0. The Yukawa potential 
tail in Eq. (25) will bring the required long-range attraction. It would be interesting to find 
out the amount of middle-range attraction brought in by two correlated or uncorrelated pion 
exchanges. 

VI. APPENDIX 

Ref. @ presented results obtained from the diagonalization in a 5-dimensional basis. For 
comparison, here we need to remove the 5th basis vector which does not have a corresponding 
one in the molecular basis. The results of the diagonalization in a 4-dimensional basis are 
given in Tables V and VI for IS = (01) and (10) respectively. 

Acknowledgements. We are very grateful to L. Wilets and L. Glozman for several useful 
comments. 
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TABLES 

TABLE I. Results of the diagonalization of the Hamiltonian (22)-(26) for IS = (01). Column 
1 - basis states, column 2 - diagonal matrix elements (GeV), column 3 - eigenvalues (GeV) in 
increasing order, column 4 - lowest state amplitudes of components given in column 1. The results 
correspond to (3 = 0.437 fm . The diagonal matrix elements Ha and the eigenvalues are relative 
to 2 m N = 1939 MeV (see text) 



State 


Ha - 2 mjv 


Eigenvalues - 2 uin 


Lowest state amplitudes 


|33[6lrJ33lp<? > 


2.616 


0.718 


-0.04571 


|33[42]o[3% s > 


3.778 


1.667 


0.02479 


|33[42] [51] FS > 


1.615 


1.784 


-0.31762 


|33[42]o[411]fs > 


2.797 


2.309 


0.04274 


|42+[6]o[3% 5 > 


3.062 


2.742 


-0.07988 


|42+[42] [33] fs > 


2.433 


2.784 


0.12930 


|42+[42] [51] FS > 


0.850 


3.500 


-0.93336 


|42+[42] [411]fs > 


3.665 


3.752 


0.00145 


|51+[6] [33] F5 > 


2.910 


4.470 


-0.01789 
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TABLE II. Same as Table I but for IS = (10) 



State 


Ha - 2 rriN 


Eigenvalues - 2 mjy 


Lowest state amplitudes 


|33[6]o[33] F5 > 


3.300 


1.083 


-0.02976 


|33[42]o[3% s > 


4.367 


2.252 


0.01846 


|33[42]o[5L> s > 


2.278 


2.279 


-0.20460 


|33[42] [411] FS > 


3.191 


2.945 


-0.04729 


|42+[6]o[3%5 > 


3.655 


3.198 


-0.07215 


|42+[42] [33] FS > 


2.796 


3.317 


0.13207 


|42+[42]o[51] FS > 


1.167 


4.058 


-0.96531 


|42+[42]o[411] FS > 


4.405 


4.459 


-0.00081 


|51+[6]o[33] FS > 


3.501 


5.070 


-0.01416 



TABLE III. Parts of the energy expectation values (GeV) of the dominant 6q state in the 
cluster model and the molecular orbital basis for IS = (01) 



Energy 


Cluster model 


Molecular orbital 




|aV[42]o[51]FS> 


|42+[42] [51] FS ) 


KE 


2.840 


3.139 


Vconf 


0.385 


0.364 


Vx 


-2.384 


-2.754 


E 


0.841 


0.749 
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TABLE IV. Same as Table III but for IS = (10) 



Energy 


Cluster model 


Molecular orbital 




|sV[42]o[51]fs> 


|42+[42] [51] F s> 


KE 


2.840 


3.139 


Vconf 


0.385 


0.364 


Vx 


-1.840 


-2.437 


E 


1.385 


1.066 



TABLE V. Results of the diagonalization of the Hamiltonian (22)-(26) for IS = (01). Column 
1 - basis states, column 2 - diagonal matrix elements (GeV), column 3 - eigenvalues (GeV) in 
increasing order for a 4 x 4 matrix , column 4 - components of the lowest state. The results 
correspond to (3 = 0.437 fm . The diagonal matrix elements and the eigenvalues are relative to 2 
m N = 1939 MeV 



State 


Diag. elem - 2 mjy 


Eigenvalues - 2 mjy 


Lowest state amplitudes 


|s 6 [6]o[3% 5 > 


2.346 


0.915 


-0.10686 


| S V[42]o[3%5 > 


2.824 


1.922 


0.08922 


|sV[42] [51] F5 > 


0.942 


2.956 


-0.98854 


|sV[42]o[4H]f5 > 


2.949 


3.268 


0.05843 



TABLE VI. Same as Table V but for IS = (10) 



State 


Diag. elem - 2 


Eigenvalues - 2 tun 


Lowest state amplitudes 


N 6 [6]o[33] FS > 


2.990 


1.453 


-0.10331 


|aV[42]o[33] F s > 


3.326 


2.436 


0.09371 


|sV[42]o[51] F5 > 


1.486 


3.557 


-0.98723 


|sV[42]o[411]f5 > 


3.543 


3.899 


-0.07694 
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